Excitation spectra of strongly correlated lattice bosons and polaritons 
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Spectral properties of the Bose-Hubbard model and a recently proposed coupled-cavity model are studied by 
means of quantum Monte Carlo simulations in one dimension. Both models exhibit a quantum phase transi- 
tion from a Mott insulator to a superfluid phase. The dynamic structure factor S{k,Lj) and the single-particle 
spectrum A{k, lo) are calculated, focusing on the parameter region around the phase transition from the Mott 
insulator with density one to the superfluid phase, where correlations are important. The strongly interacting 
nature of the superfluid phase manifests itself in terms of additional gapped modes in the spectra. Comparison is 
made to recent analytical work on the Bose-Hubbard model. Despite some subtle differences due to the polari- 
tonic particles in the cavity model, the gross features are found to be very similar to the Bose-Hubbard case. For 
the polariton model, emergent particle-hole symmetry near the Mott lobe tip is demonstrated, and temperature 
and detuning effects are analyzed. A scaling analysis for the generic transition suggests mean field exponents, 
in accordance with field theory results. 
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I. INTRODUCTION 

The Mott insulator (MI) to superfluid (SF)quantum phase 
transition in the generic Bose-Hubbard modePhas attracted a 
lot of attention in recent years due to the progress in experi- 
ments on cold atomic gases in optical lattices. 12 More recently, 
there have also been significant advances in the coherent cou- 
pling of single atoms and cold atomic gases to cavity radia- 
tion (cavity quantum electrodynamics).'^ A clean realization 
of the Jaynes-Cummings Hamiltonian has been achieved by 
coupling a superconducting qubit to a microwave cavity.!^ On 
the theory side, multi-component Bose gases coupled to light 
have e.g. been shown to support a superradiant Mott insulator 
phase with polariton condensation.^ 

In parallel, several theoretical proposals have shown the 
possibility of having a state of strongly coiTelated photons 
or polaritons in solid-state systems of coupled cavity arrays 
(also referred to as polariton models or Jaynes-Cummings- 
Hubbard models),-^ and a review of work along these lines 
has been given.l^lxhe possibiUty of preparing a system of pho- 
tons in a Mott state with one photon per site is a promising 
starting point for quantum information processing. An impor- 
tant feature shared with cold atomic gases coupled to light is 
the composite nature of the polaritons. Particularly attractive 
properties of cavity arrays would include accessibility of lo- 
cal properties in measurements and scalability in terms of size. 
Perhaps the most likely candidate for setting up such a model 
experimentally is based on extending the work on supercon- 
ducting qubits to aiTays.^ '° In contrast to cold atomic gases, 
where the interaction and/or hopping strength can be varied, 
the phase transition may be observed by changing the detun- 
ing between the two-level system and the resonator. Analy- 
sis of coupled cavity models is fruitful in its own right, as a 
detailed understanding of the corresponding models offers in- 
sight into strongly correlated polariton systems. An important 
aspect of such studies is the extent to which such systems re- 
semble the familiar Bose-Hubbard physics. 

From the above examples and many more in the litera- 



ture, it is apparent that interacting boson systems on a lattice 
are of great interest for the progress of both theory and ex- 
periment. Compared to Bose fluids, the lattice changes the 
physics in several aspects. Although long-range phase coher- 
ence still gives rise to phonon excitations — despite the break- 
ing of translational symmetry — the quenching of the kinetic 
energy makes the system much more strongly correlated.!^ 
Besides, the lattice allows the formation of incompressible MI 
states with the same integer particle number at each site. 

A large amount of work has been devoted to detailed 
studies of the Bose-Hubbard model, leading to a wealth of 
knowledge with and without additional complications such 
as trapping potentials or disorder. However, the dynami- 
cal properties and excitations in particular of the SF phase 
in the vicinity of the quantum phase transition, are still 
not completely understood. A number of authors have ad- 
dressed the dynamics of the Bose-Hub bard model in different 
dimensions, M^I'S | ^^ l " l "^ l ^'' l ^" l ^' B2Bl ,,ith results providing 

valuable information about the underlying physics, while cor- 
responding work on coupled cavity models has just begun.lSISl 
The two most important dynamic observables are the dy- 
namic structure factor and the single-particle spectral func- 
tion, which are also at the heart of theoretical and experimen- 
tal works on Bose fluids.^'' Experimentally, the dynamic struc- 
ture factor may be measured by Bragg spectroscopy or lattice 
modulation (in cold atomic gases) as well as by neutron scat- 
tering (in liquid hehum), and single-particle excitations of op- 
tical solid-state systems are accessible by means of photolu- 
minescence measurements. 

Whereas the standard Bose-Hubbard model only supports 
MI and SF phases, the physics of the polariton models is 
slightly richer. Owing to the composite nature of the con- 
served particles (polaritons), these phas es can eith er be of po- 
laritonic, excitonic or photonic characteJ ^^ ' ^^ ' ^^ ' ^^ lwith distinct 
dynamic properties. Which of the cases is realized depends 
on the value of the detuning between the cavity mode and the 
transition frequency of the atoms that mediate polariton re- 
pulsion. Very recently it has been proposed that the fractional 
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quantum Hall effect may also be realized in coupled cavity 
arrays P2l 

In general, accurate and unbiased results are very hard 
to obtain. Most existing work on spectral properties in the 
Mott phase is based on mean-field and/or strong-coupling ap- 
proximations, in which fluctuations of the particle numbers 
are more or less restricted. Results of extensive strong cou- 
pling expansions for the phase diagram'-' ^' do, however, agree 
very well with precise density-matrix renormalizati on gr oup 
(DMRG)'"* and quantum Monte Carlo (QMC) results.l?fi'?^'Bo- 
goliubov type descriptions have been found to accurately de- 
scribe the SF phase only in the limit of weak interaction, and 
fail to account for the transition to a MI and correlation fea- 
tures in the SF close to the transition. Hence the most interest- 
ing (and most difficult) regime is that near the quantum phase 
transition, where quantum fluctuations and correlation effects 
cannot be neglected. 

In one dimension (ID), quantum fluctuation effects are par- 
ticularly pronounced and mean-field methods are in general 
insufficient. Notable exceptions include situations where cou- 
pling to additional degrees of freedom provides an effective 
long-range interaction.* An interesting aspect of ID is that for 
strong (repulsive) interaction, fermions and bosons behave in 
a very similar way, and that the low-energy, long-wavelength 
physics is described by the Luttinger liquid model. 

In the present paper we employ the directed loop quantum 
Monte Carlo method,^^ which is exact and therefore yields un- 
biased results also in difficult parameter regimes. Importantly, 
our simulations preserve the full quantum dynamics. 

Few nonperturbative results are available for the spec- 
tra in the Bose-Hubbard model, namely for the dyna mical 
conductivity,'"* for the dynamic structure factor S{k, ujj^^^ 
on small systems, and for the single-boson spectral function 
A{k, (jj) in the Mott phase deduced from small systems,'^ all 
in ID. For the polariton model considered here, only A{k, oj) 
in the Mott phase has been calculated.l^S 

The focus of our work is therefore on the calculation of 
excitation spectra for both the Bose-Hubbard model and the 
polariton model within and around the first Mott lobe (i.e., the 
lobe with density one), for which comparison to recent ana- 
lytical and numerical results is made. Other issues addressed 
include the sound velocity in the SF phase, particle and hole 
masses, as well as temperature and detuning effects for the 
case of the polariton model. 

Our simulations are performed at low but finite temper- 
atures. On one hand, this complicates the analysis of the 
results, but o n the other hand it matches the experimental 
situation.lSEl 

The paper is organized as follows. In Sec. |ll]we introduce 
the two models considered. Section |lll] contains some details 
about the method. Results are discussed in Sec. |IV| and in 
Sec.|V]we present our conclusions. 



microcavities, each of which contains a single two-level atom 
with states ||) separated by energy e. Within the rotat- 
ing wave approximation one such cavity is represented by the 
Jaynes-Cummings Hamiltoniar|22l(7i — 1) 



Ht = e|T.)(T.|+^oaIa. 



+gi\U) a.|«. + U.>(T.|aI) 



(1) 



Here ujo is the cavity photon energy, and A ~ e — ujq defines 
the detuning. The atom-photon coupling g (aj, are pho- 
ton creation and annihilation operators) gives rise to forma- 
tion of polaritons (combined atom-photon or exciton-photon 
excitations). Allowing for nearest-neighbor photon hopping 
between cavities with amplitude t leads to the lattice Hamil- 
tonian 
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(2) 



ip,i, with 7ip_i ^ 



The conserved polariton number Np — J^i 
a]oL + Iti) is determined by the chemical potential 
/LtpS'Polaritons experience an effective repulsion C/eff(np) [see 



Eq. (Ill] due to the nonlinear dependence of the single-site 



energy on the local occupation number rip. We use g as the 
unit of energy and set cuo/g, fca and the lattice constant equal 
to unity. The rotating wave approximation becomes unjusti- 
fied for g comparable to e. The motivation for setting g e is 
direct comparison to previous work. The Hamiltonian (j2]i has 
been studied in II7I9II0I24I25I27I28I29I35I38I39I40L 
We also consider the Bose-Hubbard Hamiltonian 



H 



BHM 



^)~^lN, (3) 



describing soft-core bosons with repulsion U and hopping t. 
Here N — f^i — Si ^I^i' '^^e total number of bosons, 
and we use U as the unit of energy. 

As an alternative to the spin language used here, the po- 
lariton model (pj can be written as a two-band Bose-Hubbard 
Hamiltonian;""^ one boson species is itinerant, whereas the 
other is immobile (corresponding to localized excitons) with 
a hard-core constraint. This correspondence provides a direct 
connection to recent work on cold atomic gases in optical lat- 
tices, with the natural extension to the case where the excitons 
are mobile as well.* 

We shall see below that owing to the composite nature of 
the bosonic particles in the polariton model, it is generally eas- 
ier to understand the features of the Bose-Hubbard model first, 
and then explore similarities to the polariton model. More- 
over, analytical approximations are more readily available for 
the Bose-Hubbard model and provide insight into the numer- 
ical data. Periodic boundary conditions in real space are ap- 
plied in all simulations, and the system size is denoted as L. 



II. MODELS 

The polariton model we consider is the simplest among sev- 
eral recent proposals.^ ^'^^^ It describes an array of L optical 



III. METHOD 



We use the directed loop method,SI a generalization of the 
loop algorithm,"*^ '^^ which has no systematic errors and is ef- 
ficient (low autocorrelations), facilitating the simulation of 
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large systems at low temperatures. We make use of the ALPS 
library ''^ and of the ALPS applications/'' which use the 
stochastic series expansion (SSE) representatiorl^ of world- 
line path integrals. We have verified that we obtain the correct 
phase boundary in ID for selected points in parameter space. 

In contrast to most previous QMC calculations of the Bose- 
Hubbard model, the focus of the current paper is on dynami- 
cal properties. The SSE representation has the drawback that 
dynamical correlation functions in imaginary time, which we 
need to obtain spectra, are very inefficient to calculate, since 
they involve a convolution of Green functions at different SSE 
distances.'^ On the other hand. Green functions can be mea- 
sured easily in an imaginary time representation. For this rea- 
son we revert a mapping from continuous time to SSE^^when 
measuring Green functions. To each operator in a given SSE 
operator string we associate a time re [0,(3] which is stochas- 
tically sampled out of a uniform distribution. This maps the 
SSE configuration into a worldline configuration in continu- 
ous imaginary time.'^^l 

Correlation functions of diagonal operators can then be 
measured directly. For example, in the case of {pi{T)pj{0)) 
we evaluate the density pi{T) on a fine time grid. This time 
discretization limits the high energy range of the Green func- 
tion, but does not introduce any discretization error to the 
QMC algorithm itself. With the Fourier transformation of the 
density JF(pj(r)) = Pk.uj^ we measure the correlation func- 
tion J^((/5i(t)pj(0))) = {pk.,u>P-k,u) using fast Fourier trans- 
forms. 

The evaluation of off-diagonal single-particle correlation 
functions of the form (r)V'J(O)) requires some care. We 
again make use of the worldline picture, in which two opera- 
tors if)^ and ifj are inserted whenever a new loop update starts. 
Let us assume that -ip moves around (loop head) while if}^ is 
pinned (loop tail). The time and position of the loop tail are 
set as the new origin of our coordinate system, and we store 
the values (aj^/^j (r)?/'J(0)|/3) whenever the loop head ^/'^(t) 
crosses a point on the time grid with distance (i, t) from the 
new origin. Here |/3) are the states in the world line 
configuration prior to the arrival of the loop head. We then 
again use fast Fourier transformation to evaluate the correla- 
tion functions in Fourier space. 

Let us now define the observables of interest. The quantum 
phase transition can be detected by calculating the superfluid 
density p^, measured in the simulations in terms of the spa- 
tial winding number w as ps = L{'uP') / (3^^^ P — 1/kT 
being the inverse temperature. Another important observ- 
able in the context of the Ml-SF transition is the total density, 
n — {N)/L in the Bose-Hubbard model, and rip = {Np)/L 
in the polariton model. 

Concerning dynamical properties, we compute the dynamic 
structure factor S{k,uj) and the single-particle spectral func- 
tion A{k,LLj). The dynamic structure factor at momentum k 



and energy iv is given by 
1 



S{k,Lu) = 



2ttL 



— oo 



dre-- (pfe(r)pI(O)) (4) 



with the grand-canonical partition function Z and the energy 
of the nth eigenstate En. In our simulations, S{k,uj) is ob- 
tained from 



(P/c(T)/5-fe(0)> = / dujS{k,uj) 



1 + e-'^P 



(5) 



by means of the maximum entropy method. 

For the Bose-Hubbard model, the density operator pi — 
fii, and pj, — J2q ^\+k^q- For the polariton model, we can 
calculate the dynamic structure factor for photons [5p*'(/c, uj)], 
atoms [S^\k,uj)] orpolaritons [S{k,u!)] by using 



for photons , 

|Ti)(Ti| for atoms, (6) 

ajflj + Iti) (til for polaritons , 

respectively. 

The single-particle spectral function is defined as 

Mk,c^) = --Im((V'fe;^l.)). (7) 

TT 



E 



(mlV-iln) 6[uj - {E,n - En)] . 



where the real-space operator ipi entering the Green function 
is given by tpi — hi for the Bose-Hubbard model, and by V^.; = 
ai for the polariton model, maximum entropy is again used to 
map to real frequencies. 

The QMC algorithm samples the partition function in the 
grand canonical ensemble. However, using only those con- 
figurations which have a given number of polaritons enables 
us to measure observables in the canonical ensemble as well. 
Here this simple but powerful trick permits us to study the 
fixed-density phase transition which occurs in the polariton 
model as a function oft/g. 

The SSE representation requires to set a maximum boson 
number per site. In the Bose-Hubbard model, we allow a max- 
imum of six bosons per site. In the polariton model we allow 
from six (Mott insulator, fixed-density transition) up to 16 (SF 
phase) photons per site. Convergence has been monitored by 
plotting histograms of the photon number distribution, and the 
cut-offs have been chosen generously such that there was no 
truncation error. 



IV. RESULTS 

We begin with a review of the properties of the Bose- 
Hubbard model and the polariton model as they emerge 
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from previous work. Whereas a substantial literature ex- 
ists for the Bose-Hubbard model, work on the polar iton 
model began only recently, based on mean- field theory,l2Eni 
exact diagonalization,^^ the DMRG,^' the variational cluster 
approach,^'*, QMC^^^ and strong coupling theory. Our dis- 
cussion focuses on ID, and follows Fisher et al.l^and Kiihner 
etal.14'' 

The Bose-Hubbard model describes the competition of ki- 
netic energy and local, repulsive interaction. Depending on 
the ratio t/U and the density of bosons n (the system is su- 
perfluid for any i > if ri is not integer), the Bose-Hubbard 
model at temperature T = is either in a MI state or in a SF 
state. The MI is characterized by an integer particle density, 
phase fluctuations and a gap in the single-particle excitation 
spectrum. In the SF phase, we have significant density fluctu- 
ations, phase coherence, and nonzero superfluid density ps, as 
well as gapless (phonon) excitations with linear dispersion at 
small k. 

For the case of one dimension considered here, a precise 
zero-temperature phase diagram in the fi/U,t/U plane has 
been determined by Kiihner et al.,'^and these data are shown 
in Fig.[TJa). There exists a Mott lobe inside which the density 
n = 1 (higher lobes with integer n > 1 are not shown), and 
which is surrounded by the SF phase. 

There are two qualitatively different ways to make a tran- 
sition from the MI to the SF' The generic MI-SF transition 
is driven by addition or subtraction of small numbers of parti- 
cles to the incompressible MI phase, the total energy cost for 
which is given by the distance in /^-direction from the near- 
est phase boundary. Since additional particles or holes (which 
Bose condense at T = 0) can move freely, the gain in ki- 
netic energy can outweigh the interaction energy, leading to 
the MI-SF transition. Across the generic transition, which is 
mean-field like in character, the density varies continuously 
and the single-particle gap closes linearly as a function of the 
distance from the phase boundary, E„ oc S, where 6 — t — tc 
or /i — /ic is the distance from the phase boundary.^^ 

There also exists a MI-SF transition at fixed density, driven 
by the onset of boson hopping due to the increase of the ra- 
tio t/U, i.e. by quantum fluctuations. It has been shown that 
this transition occurs at the tip of the Mott lobe, and that it 
has a different universality class than the generic transition.E' 
In d dimensions, the universality class is that of the {d + 1) 
dimensional XY model, so that in ID there is a Kosterlitz- 
Thouless phase transition at the multicritical point. For this 
case, the Mott gap Eg cx exp(-const/v^t7^^) closes expo- 
nentially (i.e., very slowly) as a function of the distance from 
the lobe tip,'^ and strong deviations from the parabolic lobes 
predicted by mean- field theorjC] are observed in both strong- 
coupling''' ^' and DMRG results. Another remarkable as- 
pect of the ID case is the occurrence of multiple MI-SF tran- 
sitions along lines of constant chemical potential over an ex- 
tended range of ^ 0.2 [see Fig.[IJa)].''*^^ 

The polariton model also shows a series of Mott lobes, in 
which the polariton density np is pinned to an integer (see 
Fig.fTJb) for the phase boundaries of the Up = 1 lobe obtained 
by iJmRG'''). Even for pinned the photon and exciton den- 
sities can fluctuate. Deep in the Mott phase and for > 1, 




t/g 

FIG. 1: Zero-temperature phase diagram for (a) the Bose-Hubbard 
model and (b) the polariton model in ID. We only show the Mott 
lobes with density one. These DMRG results were obtained by (a) 
Kuhner et al.li^and (b) Rossini et al.E^ 



we can approximate the ground state by a product over sin- 
gle sites, each of which is described by the Jaynes-Cummings 
eigenstates (see, e.g., [25]) 

|7ip,-) = cos6'(7ip) |rtp, j) -sin6'(np) |np - l,t) , 
|7ip,+) = sin6'(np) |np,i) +cos6'(np) |np - l,t) , (8) 

where tan6'(np) 2g^/np/[2x{np) - A], x("-p) = 

^ g^rip + A2/4, and with eigenvalues E^{n-p) = — (/i — 
uJo)np + A/2 ± x("-p)- Hence for fixed polariton number rip, 
the ground state |np,— ) is a coherent superposition of two 
states which differ by the state of the atom (or spin) as well 
as the number of photons; this hybridization provides the con- 
nection to exciton polaritons. 

The extent of the lobes in both the /i and t directions dimin- 
ishes quickly with increasing Tip due to the reduced polariton- 
polariton repulsion Ucs{np); the t = vertical width of the 
lobes in the Bose-Hubbard model is always U. At large values 
C,t > ijj — jjL {C, being the coordination number), beyond those 
considered in the present work, the polariton model shows an 
instability.'"'' 

In this work we restrict our discussion to the region in 
the phase diagram in or close to the Mott lobes with density 
7ip — \ or n — 1. This lobe is the largest in the polariton 
model with zero detuning, and quantum effects are most pro- 
nounced. A density of on e is a lso the most interesting case for 
experimental realizations. 121211 
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All the discussion so far has been for T = 0. Both exper- 
iments and our simulations are carried out at low but finite 
temperatures, with several important consequences. Strictly 
speaking, there is no true MI at T > due to thermal ex- 
citations. However, there exist quasi-MI regions which have 
finite but very small compressibility (see also the discussion 
of temperature effects later). As long as the density remains 
close to an integer, these regions may be regarded as Mott 
insulating. Corresponding "phase diagrams" at finite T have 
been obtained for both the polariton and the Bose-Hubbard 
model.^^^"*^^ Except for our analysis of temperature effects 
in Sec. |IV] the simulations have been carried out at values of 
(3 — 3L, large enough to ensure that we have an (almost) in- 
teger density in the Mott phase. 

The Bose-Hubbard model in more than one dimension (and 
most likely the polariton model as well) exhibits a phase tran- 
sition from a SF to a normal state (gapless with no phase co- 
herence), related to the well-known A transition in liquid he- 
lium, at a temperature Tx.^*'^' This gives rise to an interven- 
ing normal region in the phase diagram, between the MI (at 
small t/U) and the SF (at large t/U)?'' In the ID case con- 
sidered here we have T\ — 0, so that for any T ^ only 
quasi-MI and normal states exist in the thermodynamic limit. 
However, when the temperature is so low that the SF correla- 
tion length in the thermodynamic limit far exceeds the system 
size L, results will be representative of the SF state. Mak- 
ing use of finite size and finite temperature effects, a scaling 
analysis in fact yields accurate results for the T = phase 
boundaries.'^^'^^'Remarkably, interacting ID bosons can be re- 
alized using cold atomic gases (the Tonks -Girardeau gas)^^E2l 
and are described by the Bose-Hubbard model at low but finite 
temperatures.-'^ 

Similar to Bose fluids, the low-energy excitations in the 
SF phase are phonons. Within Bogoliubov theory, these 
quasiparticles are described by a creation operator = 
Ufc&fc + and they have been observed experimentally 

in ultracold atom systems.!^ As some of our results can be 
understood in terms of Bogoliubov theory, let us state some 
key results for the Bose-Hubbard model. The coefficients of 
the coherent superpositions of particle and hole excitations are 
given h'^^ 



Ufc 



Vfe 



2 _K{k) + UqU + LUk 

2ujk 

2 _K{k) + naU - ujk 



2uJk 



with excitation energy 



LOk = ^K{k){2noU + K{k)), 
K{k) = 4tsin2(fc/2) . 



(9) 



(10) 



Here no is the condensate fraction, equal to no = /U in 

the simple Bogoliubov approach at T = O.'^For small fc « 0, 
we have a linear dispersion Uk ~ ±^/2notUk, and both \uk\ 
and I Vfe I are nonzero. For large fc w tt, the energy dispersion is 
±(-cfc^ + 2V4<2 + 2noUt) and thus free particle like. If we 
assume t ^ U, which is the parameter region where Bogoli- 
ubov theory is valid, then |ufcp w 1 and jv^p « for fc 3> 0, 
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FIG. 2: (color online) Single-boson spectral function A{k, lj) of 
the ID Bose-Hubbard model, for different hoppings t (and total den- 
sity n), corresponding to (a) the MI phase, (b) just below the MI-SF 
transition, (c) just above the transition, and (d) the SF phase. Here 
j-i/U — 0.5, L — 64 and l3U = 3L. Here and in subsequent spec- 
tra, the symbols and errorbars indicate the maxima of the peaks and 
the associated errors obtained by the maximum entropy method. As 
discussed in Sec. |IV A2l features with very small spectral weight are 
difficult to determine accurately. The solid red lines in (a) are mean 
field results.'^ The solid lines in (c) and (d) are the Bogoliubov re- 
sults, while the dashed lines are a fourth order approximation (see 
text).*^ 



i.e. only one excitation branch is populated at large momenta. 
This also holds true for the parameters studied numerically in 
this work. We further compare to the higher order approxima- 
tion proposed in Ref. (62^ For the Bose-Hubbard model, the 
latter yields the same equations for |ufcp, jv^p and w^,, but no 
is determined self-consistently, allowing for depletion effects. 

In the case of free bosons at T = 0, all particles condense in 
the same fc = state. However, finite temperature and/or in- 
teractions cause a certain fraction of these particles to occupy 
states of higher energy. Indeed, both for [7^0 (noninter- 
acting bosons) and no ^ (high temperature limit) we have 
|ufcp ~ 1, |vfep = 0. Moreover, with decreasing U or no, 
|vfep approaches zero most quickly at large k since in this case 
K{k) 3> uqU so that ujk ~ K{k), canceling the term —ujk 
in the expression for |vfcp. This will explain the temperature 
evolution of the single-particle spectrum shown in Sec.|IV| 



A. Bose-Hubbard model 



Despite the extensive literature on this model, there are few 
nonperturbative results available for the spectra, as mentioned 
in Sec. [l] Therefore, we investigate the single-boson spectral 
function A{k,uj) and the dynamic structure factor S{k,uj), 
with results shown in Figs. |2] and |4] 
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FIG. 3: (color online) Quasiparticle weights Uk and of the gap- 
less modes at t/U — 0.2. The symbols are integrated intensities 
from QMC and maximum entropy, the lines are the predictions from 
Bogoliubov-theory. The inset shows data at t/U — 0.14. Again, 
11 = 0.5, L = 64 and /3 = 3L. 



1. Single-particle spectrum 

Menotti and Trivedi reviewed previous work on the single- 
particle spectrum, and presented results from a random phase 
approximation.^^ Their main findings are as follows. For 
large t/U, a weakly interacting SF exists, and the spectrum 
consists of the usual two gapless phonon modes which ex- 
haust the sum rule for A{k, lo). Reducing t/U, two additional 
gapped modes appear at small k whose spectral weight in- 
creases upon approaching the quantum phase transition. At 
the transition, one of the phonon modes evolves into the par- 
ticle or hole mode (depending on which of the gaps £^g p, iJg h 
is smaller), whereas one of the gapped modes in the SF be- 
comes a gapped mode in the MI. Menotti and Trivedi^ argued 
that the appearance of gapped modes and the redistribution 
of spectral weight from coherent phonon modes to incoher- 
ent gapped modes indicate the strongly correlated nature of 
the SF state near the transition. Let us point out that particle 
and hole dispe rsions in the MI have b een calculated by several 
authors before,l '^ ' '^ ' '^ ' '^ ' '^ ' ^^ ' ^^ ' '^^ ' ^'^ ' ^ ^whereas the full spectral 
function of the MI (which also reveals the spectral weight and 
the width of the excitations) was only shown in 1 18 23 1. 

Our numerical results for the single-particle spectral func- 
tion A{k,uj) are shown in Fig. |2] The four different values 
of the ratio t/U cover the range in which the generic MI-SF 
transition takes place. According to Fig. [ija), for the chosen 
value of ji/U — 0.5 the transition occurs at t/U k, 0.14. 
In each panel we also report the total density n to three dec- 
imal places, although our simulations provide much higher 
accuracy. The MI [(a) and (b)] exhibits the familiar gapped 
particle and hole bands.'Sl The additional particles exhibit a 
free-particle dispersion since the energy penalty for double 
occupation is the same at every site. In particular, we see in 
Fig. |2|a), (b) that the particle band width is St (the factor of 
two arising from the fact that particle hopping involves a dou- 
bly occupied site), whereas the hole bandwidth is 4i. The 
Mott gap decreases with increasing t and a symmetry of par- 
ticle and hole bands emerges.' In addition to our QMC 



results we plot the mean-field dispersiorP^ in Fig. |2]^a). For 
larger t/U — 0.13, mean-field theory already predicts a su- 
perfluid, although the critical hopping in ID is tt./U « 0.14. 

In the SF phase [Fig. |2|c),(d)], we obtain the expected 
Goldstone modes with linear dispersion at small k. Addi- 
tionally, we see two gapped signals whi ch we relate to the 
gapped modes discussed by other authors .l^^lilllll whereas the 
negative-energy gapped mode is clearly visible in Fig. [2|c) 
just above tc, the gapped modes have almost disappeared in 
Fig. |2]^d). Since we approach the phase transition above the 
lobe tip iii/U — 0.5) the particle band becomes the gapless 
mode and carries more spectral weight, while the gapped hole 
band evolves into a gapped mode in the SF. This agrees well 
with the findings of Menotti and Trivedi .l^^l In accordance with 
Bogoliubov theory, the excitations in the SF phase are free- 
particle like for large k. The bandwidths of the excitations 
both in the MI and the SF phase scale roughly linearly with t. 

In Figs. |2|c) and (d) we also show results for the phonon 
dispersion ±ujk (without taking into account the weights |ufe|, 
|vfc|) from Bogoliubov theory as well as the higher-order ap- 
proximation of Ref. 62 Whereas the simple Bogoliubov ap- 
proach (neglecting depletion of the condensate) agrees quite 
well with our data despite the rather small value of t/U, we 
do not find the higher order approach to be systematically bet- 
ter. In particular, at large k, the phonon bandwidth is notice- 
ably underestimated, which may be a result of an overestimate 
of depletion effects (these are most visible at large fc). The 
agreement with Bogoliubov theory at small k coincides with 
the findings of Menotti and Trivedi. Rey et al.^- found the 
higher order approximation to be consistent with numerical 
results for other observables but do not show the spectra seen 
in Fig. |2] Note that these authors consider larger particle den- 
sities n > 5 where the Bogoliubov-type approximations are 
more reliable. Finally, we tried to use our QMC results for 
the superfluid fraction for no in the expressions obtained from 
Bogoliubov theory, but the results are worse than for tiq — n. 

The spectral weight of the excitations decreases with in- 
creasing k in all spectra of Fig. |2] although this is more pro- 
nounced in the SF phase than in the MI. In Fig. |3] we show 
the quasiparticle weights of the massless modes in the SF 
phase, obtained by integrating over the quasiparticles peaks in 
the spectra, and compare them to Bogoliubov theory (Eqs.[9] 
and 10 1. We verified that the QMC spectra satisfy the sum 



rule. The spectral weight of the lower branch decreases more 
quickly, consistent with the Bogoliubov picture. However, 
Bogoliubov theory overestimates the quasiparticle weights, 
especially at small fc. Besides, there is a significant broaden- 
ing of the peaks on approaching the zone boundary. At strong 
coupling close to the phase transition (inset of Fig. [3]l, the 
quasiparticle weight of the lower branch decays much more 
quickly than Bogoliubov theory would predict. 

Sound velocity. The sound velocity — ^j^lfe^o of the 
phonon excitations in the SF phase was calculated for the 
Bose-Hubbard model by Menotti and Trivedi using a random 
phase approximation.^^ They concluded that vanishes at 
the generic transition, but remains nonzero when crossing the 
multicritical point.^^^l^li jjj fjjgj]- results, there is a very sharp 
downturn of toward zero close to t^. We are not aware of 
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FIG. 4: (color online) Dynamic structure factor S{k, lj) of the Bose- 
Hubbard model for the same parameters as in Fi g. [2| Panel (d) in- 
cludes the same analytical approximations as Fig.l2td). 

any calculations of for the polariton model. 

From our QMC simulations, we can determine from lin- 
ear fits to the spectrum. Apart from the limited accuracy of 
the maximum entropy inversion, this works quite well away 
from t(,. In agreement with Bogoliubov theory, we find for the 
Bose-Hubbard model a linear dependence oc | i — tc | and 
good agreement of results for L — 32 and 64. 

Determining the behavior of as t ^ is more difficult 
for two reasons. First, the phonon spectrum becomes nonlin- 
ear due to finite-temperature effects (see discussion below), 
rendering linear fits ill-defined. Second, the position of the 
phase transition changes with system size, so that no reliable 
finite-size scaling of can be carried out. The situation is 
similar for the polariton model, and we therefore do not show 
results for here, leaving this as an interesting issue for fu- 
ture work. 



2. Dynamic structure factor 

The single-particle spectral function provides information 
about the energy and lifetime of particles or holes added to 
the interacting ground state. In contrast, the dynamic structure 
factor — corresponding to the imaginary part of the dressed 
particle-hole propagator — yields insight into the density fluc- 
tuations in the ground state. In general the two quantities 
do not exhibit the same features. However, for broken U(l) 
gauge symmetry in the SF phase, they are both dominated by 
the same single-particle excitations (phonons).^*" We find this 
statement to hold in ID even though no symmetry breaking 
occurs. 

The density operator in the Bose-Hubbard model is = 
J2i e^^^^'fii. For fc = 0, we have pj, = J^i S{k, uj) has 

a trivial contribution at = which we dismiss by consider- 
ing pI. = J2i e~''^'(ni — {fii})- The above-mentioned relation 
to particle-hole excitations becomes evident by rewriting the 



density operator as p[ = J^q Kj+k^g- 

We show results for S{k, uj) in Fig.|4] According to Huber 
et al.^ S{k, uj) in the MI phase should exhibit a continuum of 
paiticle-hole excitations, starting at uj ^ E„ due to the Mott 
gap in the single-paiticle spectrum (see Fig. [2]). For the pa- 
rameters in Fig. Ra), Eg/U ~ 0.7. The dispersion of the 
particle and hole oands is very weak. Note that we find no 
agreement with the two single-particle excitations Eg+eh{k), 
Eg + ep(fc) discussed by Huber et al. This may be a result of 
their mean-field treatment of the two-dimensional case. Our 
results do agree qualitatively with exact numerical results on 
small clusters.Q21 

For larger t/U, the Mott state contains nontrivial density 
fluctuations, and the upper band in S{k,uj) acquires some k 
dependence. The energy of the excitations in S{k, uj) [fol- 
lowing X)q{^h('?) + ep(fc ~ 9)}]'^ generally increases with in- 
creasing k. This is obvious from the momentum dependence 
of the particle and hole bands in A{k, uj), and also agrees with 
the expectation that long-wavelength density fluctuations in 
a Mott state require less energy than fluctuations with short 
peiiods in real space. 

For t ^ tc in Fig. Qb), we find a low-energy mode with 
nonlinear dispersion, which we interpret as a precursor of the 
linear excitations of the SF phase [see panel (d)]. Even for 
t }ttc [Fig.|4]^c)], the gapless low-energy mode in our numer- 
ical results is not linear A linear spectrum is a result of the 
condensation of bosons in the SF phase, but is not expected 
in the normal phase. Since our simulations are done at finite 
temperature, and because the phase coherence length is small 
close to tc, we can understand the absence of a clear, linear 
signature in Fig.Qc). Going to larger tc, we indeed see linear 
excitations near fc = [Fig.|4|d)]. Similar effects are expected 
for the single-particle excitations, but are difficult to see on the 
scale of Fig.|2] Coming back to Fig.|4]^c), away from fc = 0, 
we find a free-particle like contribution, similar to the case of 
the MI. This excitation carries negligible spectral weight near 
fc = 0. 

Apart from finite-temperature effects, these features are 
qualitatively similar to the excitations discussed by Huber et 
al.,E21namely a gapless sound mode (related to phase and den- 
sity modulations) dominant at small fc, and a massive mode 
(corresponding to exchange between condensate and noncon- 
densate at fixed density) acquiring spectral weight at fc > 0. 
Additionally, we see in Fig. |4|c) the (weak) signature of a 
gapped mode at small fc, the nature of which we cannot deter- 
mine from our present simulations. For t/U = 0.2 [Fig.Qd)] 
the excitation "band" in S{k,uj) follows closely the Bogoli- 
ubov mode, in accordance with the discussion at the beginning 
of this section. 

At this point, a comment concerning the accuracy of the 
spectra obtained from the maximum entropy inversion is in or- 
der The spectral weight of the features visible in density plots 
such as Fig.Qd) varies over orders of magnitude. Some very 
weak signals, such as the group of points located at around 
k — 7t/2 below the main excitation band (with a weight that 
is a factor 10000 smaller than that of the dominant features), 
are expected to be artifacts. We shall see below that in the 
polariton model, there actually exist real excitations with very 
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FIG. 5: Single-photon spectral function A'*'{k, lo) of the ID polari- 
ton model at n/g = 0.4 for different hoppings t, corresponding to (a) 
deep in the MI, (b) just below the MI-SF transition, (c) just above the 
transition, and (d) in the SF phase. Here /3g = 3L and L = 64. With 
increasing t, the density plots are more and more "overexposed" to 
see less dominant features. 



small spectral weight which are easy to miss in the maximum 
entropy inversion. To reliably study such excitations, analyti- 
cal approaches (if available) are clearly superior.'25l 

Our findings for the dynamic structure factor are consis- 
tent wi th pre vious numerical results on small systems {L — 
can confirm the broadening of the excitations 
with increasing k in the SF phase, related to two-particle 
continua.E3 However, the maximum entropy method is not ca- 
pable of resolving fine structures as (generically) seen in exact 
diagonahzation results for small clusters.'^ 



B. Polariton model 

For the polariton model, the only published results on dy- 
namic properties are for the singl e-particle spectrum of the 
MI phase at zero temperature.'^^'^ As pointed out before, the 
nature of the conserved particles in the polariton model is de- 
termined by the detuning. We start by discussing the case 
A = for which the polaritonic character of the excitations 
is most pronounced. This can readily be seen from Eq. (l8]l. 



where [rip, |) and 



1, t) contribute with equal weight. 



1. Single-particle spectrum 

In Fig. |5] we show our QMC results for the single-photon 
spectral function. As for the Bose-Hubbard model, the values 
of the ratio t/g range from deep in the Mott phase across the 
generic transition well into the SF phase. According to a fi- 
nite size scaling analysis for /i/g = 0.4, the phase transition 
occurs at t^/g = 0.0626(1) (see Fig. 14 1, in agreement with 
Fig. [ijb). Hence panels (a) and (b) are for the MI regime. 



FIG. 6: Single-photon spectral function in the Mott phase for the 
same parameters as in Fig.|5] showing additional excitations at higher 
energies. 
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FIG. 7: (color online) Quasiparticle weights and Vfc of the gapless 
modes of the polariton model, similar to Fig.|3] 



whereas (c) and (d) are for the SF phase. 

The results in the MI shown in Fig. |5ja), (b) agree well 
with previous numerical work.'^i' Similar to the Bose-Hubbard 
model, there exist particle and hole bands, separated by the 
Mott gap. It is important to stress that although we add bare 
photons to the system, the particle and hole excitations re- 
flect the properties of the polaritons in the system. Whereas 
the ratio of particle and hole bandwidths is two to one in the 
Bose-Hubbard model, it depends on the character of the quasi- 
paiticles (polaritons) in the polariton model and varies with 
detuning. With increasing t/g, the gap closes and the band- 
widths of excitations increase (effective masses decrease). 

Recent analytical work revealed the existence of so-called 
upper polariton modes at higher energies, which represent an 
important difference between the Bose-Hubbard model and 
the polariton model. For the Mott lobe with Up — 1, only one 
such (particle) band exists, corresponding (for small enough 
t/g) to a transition between the ground state |np = 1, — ) and 
the state jrtp = 2, +) (see Eq. (jsj). The weight of this high- 
energy excitation is very small compared to the dominant par- 
ticle and hole modes discussed above (0.04 as compared to 
1.46 for the fc = atomic-limit results in f25 |); with increas- 
ing t/g the weight difference becomes even larger The en- 
ergy splitting between the — and + branches of eigenstates 
increases further for detuning A (Fig. 2 in |7 1). 

The upper polariton mode is not visible in Figs.jSja) or (b). 
Excitations with small spectral weight are notoriously difficult 
to see using QMC in combination with maximum entropy. In 
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FIG. 8: Single-photon spectrum A^^{k,u]) of the polariton model 
along the line = 1 crossing the Kosterlitz-Thouless transition. 
Here i = 64 and Bq = 3L. 



the present case, this is aggravated by the fact that the reso- 
lution of maximum entropy decreases at high energy. Nev- 
ertheless, we see a signature of the upper polariton band in 
Fig. |6] and the latter is also present (but not shown) in the 
high-temperature data of Fig. 12 a); high-energy features are 
easier to resolve in QMC/maximum entropy at higher temper- 
atures. From the eigenvalues of the states (|8]l we can deter- 
mine the excitation energy of the upper mode in the atomic 
hmit as /g = —{^ — uJo)/g + (\/2 + 1) « 3 for the pa- 
rameters of Fig.|6j in reasonable agreement with our results in 
Fig.[6|a) given the ill-conditioned nature of the problem under 
consideration. Note that the upper polariton mode can be seen 
even close to the phase transition in Fig. |6|b). The weight 
of the upper mode in Fig. |6]is about a factor of 100 smaller 
than that of the conventional particle and hole excitations. Al- 
though the upper polariton mode exists also in other results 
for the single-particle spectrum in the Mott phase (Figs. [8 12 
and[T3]l, we focus on the low-energy conventional modes with 
large spectral weight. The latter can be determined accurately 
from our simulations, and will be the dominant feature in ex- 
periments. 

Figures |5|c) and (d) contain the first spectra of the polariton 
model in the SF phase. There is a clear signature of the gapless 
phonon modes starting at A: = w = 0, with linear dispersion 
at small k. In the SF phase but close to the transition, we see 
an additional gapped mode at w < [Fig. [5]^c)]. Our results 
at these and at further couplings t/g (and t/U) suggest that 
these gapped modes disappear more quickly with increasing 
t/g than for the Bose-Hubbard model, which can be explained 
in terms of the photonic SF expected for the present parame- 
ters (see below).'^Note that a simple Bogoliubov type theory 
for the polariton model does not exist, due to the composite 
nature of polariton excitations. 

Figure |7] shows the quasiparticle weights. The general 
shapes resemble the Bose-Hubbard model case (Fig. |3]l, but 
the lower branch decays very quickly in the polariton model 
even at t/g = 0.15 quite far from the phase transition. This 
may be attributed to the fact that the energy cost for particle 
and hole excitations is different due to the dependence of [/gff 
on rip. Again, most of the spectral weight in the SF phase is 
found at small k. 

Fixed density. In Fig. [8] we show the single-photon spec- 
trum across the fixed density transition (rip = 1), obtained 
by selecting configurations only at that density. In Ref. |40, 
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FIG. 9: (color online) Effective particle (m^) and hole masses (m~) 
along the line rip = 1, as obtained from fits to the bands in AF'^{k, uj) 
near fc = 0. 



(a) I 



|]]linini]!iii! 




t/g = 0.01 
n„= 1.00 





(b) 






A 
















t/g = 0.061 




n,, = 1 .02 




t/g = 0.15 
n„ = 2.93 



FIG. 10: Polariton dynamic structure factor S{k,uj) for the same 
parameters as in Fig.|5] The insets show an extrapolation of the Mott 
gap at small fc to L — > oo. 



the critical hopping was determined as tc/g — 0.198 [cf 
Fig. [TJb)]. The spectra in both the MI and the SF look very 
similar to those across the generic transition shown above. 
This may be different very close to the multicritical point, but 
this regime is most demanding numerically if the results are 
to be used in a maximum entropy inversion. 

From the spectra obtained at constant density in the MI, we 
can estimate the effective particle and hole masses by fitting a 
quadratic dispersion to the bands in the vicinity of A: = 0. In 
the Bose-Hubbard model, there is an emergent particle-hole 
symmetry on approaching the lobe tip,GESI and similar behav- 
ior is suggested by the evolution of the particle and hole bands 
with increasing t/g also in the polariton model. For fixed po- 
lariton density, the two masses approach each other and van- 
ish at the phase transition. This has been demonstrated in 2D 
based on a strong-coupling approach. In the region not too 
close to the phase transition, where stable fits can be obtained. 
Fig. |9] confirms this observation also in ID. 
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2. Dynamic structure factor 

The evolution of the polariton dynamic structure factor 
S{k, Lu) across the MI-SF transition is shown in Fig. 10 Re- 
markably, the results look very similar to those for the Bose- 
Hubbard model. Close to the atomic limit [t/g ~ 0.01 in 
Fig. [TO|a)] we see a gapped, almost flat feature with energy 
u! « 0.6g. A look at the corresponding single-particle spec- 
trum in Fig.|8ja) reveals that this value is identical to the Mott 
gap. The almost flat particle and hole bands cause a very 
weak dispersion also for the particle-hole excitations visible 
in S{k,uj). It is useful to remember that it is the effective 
polariton-polariton repulsion mediated by the atom-photon 
coupling that determines the Mott gap. For a single site and 
Up — 1 (i.e., for the case of adding a second polariton). 



t/eff(l) = 2^g^ + (A/2)2- + (A/2)2-A/2. (11) 

For zero detuning (A = 0), Ucs{l)/g = 2 - ^2 ~ 0.59. 

As for the Bose-Hubbard model, the excitations in S{k, us) 
acquire a noticeable dispersion with increasing t/g, and the 
k = gap closes. Figures [TO|b) and (c) are both close to the 
phase transition. An inspection of the k = region shows a 
weak linear mode with very small slope 0(0.01), correspond- 
ing to the small superfluid density existing in both L = 64 
systems. The massive mode extends to fc = with a tiny in- 
tensity (thus not visible in the figure). An extrapolation of the 
gap to L = oo (insets) shows that it scales to zero in Fig.[TO|c), 
but stays finite at the smaller hopping in Fig. 10 b). Indeed, 
a finite size scaling of the superfluid density (discussed later) 
implies that Fig. [TO{b) is just below the phase transition. In 
addition to finite size effects, we again observe finite temper- 
ature effects in the form of deviations from the expected lin- 
ear spectrum close to (see discussion for the Bose-Hubbard 
model). For even larger t/g, the spectrum exhibits a single lin- 
ear mode at small k. Similar to the spectral function, gapped 
modes seem to be suppressed quickly in the SF phase. 

Clearly, the polariton dynamic structure factor S{k, uj) rep- 
resents a useful probe to distinguish between the Ml and the 
SF phases. We have argued before that the polariton Ml has 
fluctuations in the photon and exciton density, whereas the po- 
lariton density is pinned. 

We now demonstrate that the exciton (atom) and photon 
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structure factors, S'^^'{k,uj) and S^^{k,LLj), shown in Fig 
do not reflect this fact, and therefore cannot be used to char- 
acterize the nature of the Mott state. To this end, it is im- 
portant to notice that the Jaynes-Cummings Hamiltonian has 
two branches of eigenstates \np, +) and |np, — ) (the latter con- 
taining the ground state, see also Eq. ([8])) with the same po- 
lariton number but different energy.^^ In the atomic limit, the 
dynamic structure factor for photons and excitons [Eq. (|5]l] is 
dominated by the contributions (rip, — | pj, |?^p, — ) (with uj — 
0) and (rip, +| joj. Jnp, — ) (with energy uj — 2^ng'^ + A'^/A, 
equal to 2 in Fig.[Tl|. Any additional peaks at finite t/g have 
much smaller spectral weight and cannot be accurately re- 
solved by our method. However, since the matrix elements 
of exciton and photon density operators for the combination 
(np, + \ pI |«p, — ) have the same modulus but opposite sign. 
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FIG. 11: Dynamic structure factor for excitons (5"'') and photons 
j^^ph-j ^^j, jjjg same parameters as Fig. 10 a). 
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FIG. 12: Temperature dependence of A^^{k,uj) in the polariton 
model at n/g — 0.4, in the MI (left) and in the SF phase (right). 
Here again L = 64. Corresponding results at I3g = 192 can be 
found in Figs.|5|a),(d). 



the dominant contributions to S^*^{k,uj) and SP^{k,Lu) can- 
cel in the case of the polariton structure factor S{k,uj). The 
dispersionless excitations near u/g = 2 seen in Fig. [TT] are 
therefore absent in S{k, uj), as confirmed by our data. Hence, 
while the upper polariton modes in the single-particle spec- 
trum have small but finite weight, here their contribution is 
zero. Consequently, the polariton dynamic structure factor 
closely resembles S{k, uj) of the Bose-Hubbard model. 



3. Temperature effects 

Experimental reaUzations of Bose-Hubbard models using 
cold atomic gases are usually prepared very close to zero tem- 
perature (nK range). In contrast, due to the strong matter-light 
coupling achievable in cavities, reaUzations of polariton mod- 
els offer a chance of operation at significantly higher temper- 
atures. The critical temperature at which the Ml state starts 
to lose its characteristic integer density has b een estimated for 
the polariton model as T*/g « O.OS.I^^l] 

For feasible val- 
ues of the coupling g, T* falls into the mK range. Gener- 
ally, Mott-like physics is expected as long as the Mott gap is 
significantly larger than the thermal energy (and the number 
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of particle-hole excitations is small). The finite-temperature 
physics of t he Bose-H ubbard model has been analyzed by sev- 
eral groups i^l ^ l ^^ l ^^ I Here we consider the effect of low but 
finite temperatures on the excitation spectra of the polariton 
model. This also provides information about the sensitivity 
of the results to the (necessarily finite) value of j3 used in our 
simulations. 

The present method permits calculation of spectra also out- 
side the Ml, i.e., in the SF and the normal phase. We have 
pointed out above that, strictly speaking, there is no SF phase 
at T 7^ in ID. Nevertheless, SF like properties can be seen 



for T sufficiently small. The results in Fig. 12 underline the 
discussion of finite-temperature effects on the dispersion in 
the SF near k — Q. The particle excitation is obviously not lin- 
ear in panels (b) and (d), for which the temperature is higher 
than in Fig. |5] 
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At finite but 



Results for A^'^{k,Ld) are shown in Fig. 
low temperature [Fig. [T2|c,d)] they still closely resemble the 
results at T w in Figs. [5|a) and (d). At high temperature 
[Fig. [T2|a,b)], we observe strong broadening of the particle 
band at all fc, and strongly suppressed spectral weight for hole 
excitations. Existing work for the Bose-Hubbard model finds 
that at finite temperature additional multi-particle and hole 
bands arise.*''' We see an additional excitation for (3g = 4.4 
and t/g = 0.01 at an energy of oj/g ~ 3.1. The weight of 
that excitation is about 50 times smaller than the weight of 
the main peak with energy uj/g ~ 0.2 and thus not shown in 
Fig. 12 (a). This excitation is consistent with the upper polari- 
ton mode discussed above. 

We note that a broadened "gapped" spectrum is compatible 
with a density that deviates from the integer value character- 
istic of the Ml, and this has to be kept in mind for potential 
applications relying on integer density. The numerical results 
for the total density are shown in each of the panels, demon- 
strating that despite the large particle-hole gap the polariton 
density deviates significantly from the low-temperature value 



rip = 1 for the parameters of Fig. 12 a). 

Some of the features observed in the SF phase can be ex- 
plained by means of Bogoliubov theory for the Bose-Hubbard 
model. In particular, we have discussed above that with in- 
creasing temperature (where the condensate fraction uq 0) 
the spectral weight of the negative energy branch vanishes first 
at large k, in agreement with our numerical results. In ad- 
dition the broadened positive energy branch no longer has a 
clear linear behavior at small k. 

We would like to point out that not only finite temperature 
but also disorder is an inevitable feature of experimental real- 
izations of coupled cavity systems. Although not studied here 
directly, it has been stated^^ that the effect of disorder (in the 
form of local variations of the parameters uiq, g and t) has 
similar consequences as finite temperature. 



4. Detuning 

The detuning between the cavity photon mode and the 
atomic level splitting is an important parameter in the polari- 
ton model which is absent in the Bose-Hubbard model. Its 
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FIG. 13: Single-photon spectra with detuning A = e — ojo for the 
excitonic case A/g — —2, /i/g = —0.5 (a,b) and photonic case 
A/g = 2, n/g = 0.64 (c,d), in the MI (a,b) and in the SF (c,d). 
Here i = 64 and 3q = 3L. 



influence on the physics has been discussed before.Sl22153|j)g. 
tuning can also be easily changed experimentally, motivating 
a calculation of the excitation spectra for A ^ 0. Our results 



are shown in Fig. 13 



The extent of the different phases, namely excitonic or po- 
laritonic Ml and photonic or polaritonic SF in the phase dia- 
gram has been analyzed for a two-site system.^^ The way to 
distinguish between the polaritonic SF and the photonic SF 
is to monitor fluctuations in the exciton occupation number 
(pinned in the photonic SF but fluctuating in the polaritonic 
SF). The conclusion has been that for t « |A| a polaritonic 
SF exists only for A/g < —1 This would match with the con- 
jectured photonic nature of the SF near the lobe tip in two 
dimensions.^'' However, it is not clear if these strict values 
also hold for larger systems and the thermodynamic limit. Be- 
sides, the work by Irish et alW^ is exclusively concerned with 
the fixed density transition occurring at t w |A|, whereas a 
polariton SF may exist also for t < |A| if density fluctuations 
are allowed (generic transition). The spectrum in Fig.fljj^b) is 
for such a set of parameters. 

Again pertaining to the fixed density case, the Ml state is 
supposed to be of excitonic nature for A/g < — 1, t < |A|, 
and of polaritonic nature for \A\/g < 1 and small enou gh t /g 

a). 
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respectively </ 1 A | .^^ The former case is depicted in Fig 
whereas the latter corresponds to the A = results reported 
inFig.|5] 

For A ^ g, photon excitations are always lower in energy, 
and the effective interaction approaches zero. As a result. Ml 
regions are very small or nonexistent, and the photonic SF 
state is similar to that of the Bose-Hubbard model in the limit 
of large t/U.^'* 

Here we consider A/g ^ ±2 for comparison to previous 
calculations of the spectra in the Mott phase. (Note that 
the rotating wave approximation formally requires |A| <C 
e,ujQ.^'^) These correspond to effective repulsions Ucs/g = 
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0.096 (for A/g = 2) respectively Ucs/g = 2.096 (for 
A/g = —2), in excellent agreement with the width of the 
Ti„ — 1 Mott lobes for the same parameters.l^S 



Our results in Fig. 13 show that again the spectra are dom- 
inated by the generic features of the Ml and the SF. However, 
the detuning in the present case changes the ratio of the band- 
widths of particle and hole bands W^/Wh in the Mott state.l^ 
While for A = 0, W^/W^ ~ 3, we find Wp/M^h ~ 2 (sim- 
ilar to the result for the Bose-Hubbard model) for A/g = 2 
and Wp/Wh ~ 7 for A/g = —2. The incoherent features 
observed for A/g ~ —2 in Ref 24 are not seen here. As 
mentioned before, the energy of the upper polariton modes 
(not shown) increases for A OJ 

In the SF, we find the expected gapless excitations, as well 
as gapped modes indicative of a correlated superfluid. Since 
for A /g = 2, [/gff is very small, the Mott gap of the dispersive 



bands in Fig. 13 c) is also small (0.057 g), but it is still larger 
than the temperature scale in our simulation T/g = 0.005. We 
note that, within our resolution, the positive energy spectrum 
in Fig. [T3jd) looks gapless, but not clearly linear. In this re- 
spect, the spectra for finite detuning resemble those obtained 
at high temperatures. Apart from this issue and a scaling of 
energies (due to the dependence of Ues on A), the spectra ob- 
tained for A/g = —2 are very similar to those for A = 0, 
whereas those for A/g = 2 resemble closely the results for 
the Bose-Hubbard model. 



5. Phase transition 

To end with, we present a scaling analysis for the generic 
phase transition. As pointed out by Fisher et al.E the scaling 
relation 



^p{5L^'\(i/L^) 



(12) 



should hold for the superfluid density across the Ml-SF tran- 
sition. Here v is the critical exponent of the correlation length 
which is expected to diverge like ^ ^ (5^", and z is the 
dynamical critical exponent. The generic transition in the 
Bose-Hubbard m odel has mean-field exponents z = 2 and 
V ^ l/z ^ l/2.CEiIl2ll! Recent field theor>nSl and strong- 
coupling results^^ predict the same universality classes for the 
polariton model, in conflict with numerical results in two di- 
mensions which suggest the absence of multicritical points.^-' 
We test in ID the scaling hypothesis Eq. ( 12 1 with z = 2 



and the hyperscaling relation z = 1/v, along the line [i/ g ~ 
0.4 where the generic transition is expected [see Fig. [TJb)]. 

To this end we keep the temperature constant at /3 = L^/10 
and plot psL'^^^^'^ over {t — tc)L^^'^ to obtain the universal 



function p (Fig. 14 1. Defining a cost function in the spirit of 
Ref 70, allows us to evaluate the quality of the finite-size- 
scaling plot quantitatively. 

We find the minimum of our cost function at <c = 
0.0626(1) and v = 0.50(2) when matching p close to the 
phase transition [\{t - tc)L^^'^\ < 1 in Fig.[l4j. We note that 
with the system sizes available, this result is not very stable. 
A fit in a larger region of |(t — tc)L^^'^ \ provides a better data 
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FIG. 14: (color online) Finite size scaling for the generic transition 
in the polariton model at /i/g — 0.4, testing the scaling hypothe- 
sis Eq. l|T2j. 



collapse overall (but worse close to the phase transition), with 
f « 0.65 and very slightly smaller tc- 

A similar scaling with z = 1 did not succeed, so that 
we conclude that the universality class of the generic tran- 
sition in the polariton model is the same as that in the Bose- 
Hubbard model, despite the composite nature of the quasipar- 
ticles. This is consistent with recent field-theory and strong- 
coupling results."°'^^' 

An accurate scaling analysis for the fixed density transition 
through the lobe tip has been found to require much larger 
system sizes and is therefore not shown. In the ID case con- 
sidered, the shape of the lowest Mott lobe in ID (see Fig.[T]) 
suggests that the similarity to the Bose-Hubbard model holds 
also in this respect, i.e. a Kosterlitz-Thouless type phase tran- 
sition. 



V. CONCLUSIONS 

We calculated the single-boson spectral function and the 
dynamic structure factor of the Bose-Hubbard model, and for 
a recently proposed model of itinerant polaritons in coupled- 
cavity arrays. These models undergo a quantum phase transi- 
tion from a Mott insulator to a superfluid state upon increasing 
the hopping integral of the bosons respectively photons with 
respect to the interaction. Results in one dimension, within 
and close to the Mott lobe with density one, have been ob- 
tained. 

Despite the generally different nature of the conserved par- 
ticles, the models exhibit very similar spectral properties, in- 
cluding gapped particle and hole bands in the Mott insulat- 
ing phase, and Bogoliubov type excitations in the superfluid 
phase. Additional excitations related to the second branch of 
upper polariton states exist in the single-particle spectrum of 
the polariton model,'^but cancel out in the dynamic structure 
factor In general, these features have high energy and very 
small spectral weight, so that for practical purposes the ex- 
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citation spectra are qualitatively similar to the Bose-Hubbard 
model. 

Correlation effects are particularly strong in the one dimen- 
sional case considered. Our results in the superfluid phase rep- 
resent the first unbiased nonperturbative spectra for A{k,uj) 
in both models and for S{k,uj) in the polariton model (in 
both phases). Good qualitative agreement with recent analyt- 
ical work on the two-dimensional Bose-Hubbard model was 
found, and we have compared our results in the superfluid 
phase to Bogoliubov theory. The limiting cases of the Mott 
insulator close to the atomic limit, as well as the weakly in- 
teracting superfluid are described quite well by analytical ap- 
proximations, whereas in the phase transition region, our non- 
perturbative results show considerable deviation. Emerging 
particle-hole symmetry on approach of the multicritical lobe 
tip has been demonstrated for the polariton model. 

For the polariton model, we have also explored the influ- 
ence of detuning and finite temperature on the spectral prop- 
erties, and have presented a scaling analysis to determine the 
universality class of the generic phase transition. Keeping in 
mind experimental realizations of coupled cavity arrays, inter- 
esting open issues for future work include the excitation spec- 



tra in the two-dimensional case (and comparison to analytical 
results^^), the behavior of the sound velocity across the phase 
transition (also for the Bose-Hubbard model) and disorder. 

The present work further highlights the fact that the physics 
of strongly correlated bosons as described by the Bose- 
Hubbard model may be observed in terms of optical models 
that, if realized, would have some distinct experimental ad- 
vantages and further contain new degrees of freedom due to 
the mixed nature of the quasiparticles. 
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